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Smooth-throat wormholes are treated on as possessing quantum fluctuation energy with scalar 
massive field as its source. Heat kernel coefficients of the Laplace operator are calculated in back- 
ground of the arbitrary-profile throat wormhole with the help of the zeta-function approach. Two 
specific profile are considered. Some arguments are given that the wormholes may exist. It serves 
as a solution of semiclassical Einstein equations in the range of specific values of length and certain 
radius of wormhole's throat and constant of non-minimal connection. 
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^ ; I. INTRODUCTION 

u 

Great interest towards the space-time of worniholes dates back at least to 1916 i^l]. Subsequent activity was initiated 
' by both classical works of Einstein and Rosen in 1935 in the context of black hole space-time structure and the later 
series of works by Wheeler in 1955 with his excellent idea to create everything from nothing. The more recent 
• interest in the topic of wormholes has been rekindled by the works of Morris and Thorne ..12j and Morris, Thorne, 
and Yurtsever jl^l who made use of the concept of wormhole in scientific discussion of " time machine" . These authors 
constructed and investigated a class of objects they referred to as "traversable wormholes". Their work led to a flurry 
^ ' of activity in wormhole physics 16] . 

It is well-known that the central problem of the traversable wormholes connects with unavoidable violations of the 
null energy condition. It means that the matter which should be a source of this object has to possess some exotic 
properties. For this reason the traversable wormhole can not be represented as a self-consistent solution of Einstein 
■ equations with usual classical matter as a source because usual matter is sure to satisfy all energy conditions. One 
' way out is to use quantum fields in frameworks of semi-classical quantum gravity. The point is that the vacuum 
average value of energy-momentum tensor of quantum fluctuations may violate energy conditions. A self-consistent 
' wormholes in framework the semiclassical quantum gravity have been studied in Ref. [l4'|. In our recent paper jll| 
> , we have considered a possibility for self-consistent solution of semi-classical Einstein equations for specific kind of 
^' wormhole - short-throat flat-space wormhole. The model represents two identical copies of Minkowski space with 
^ spherical regions excised from each copy, and with boundaries of these regions to be identified. The space-time of 
this model is flat everywhere except a two-dimensional singular spherical surface. The vacuum average of energy 
of quantum fluctuations of massive scalar field with non-minimal connection serves as a source for this space-time. 
Due to the fact that this space-time is flat everywhere a complete set of wave modes of the massive scalar field can 
be constructed and ground state energy can be calculated. In the paper we present a calculation of full energy of 
quantum fluctuations rather then energy density and use the Einstein's equations with quantum source only, without 
classical contribution. We found that the energy of fluctuations as a function of radius of throat a may possess a 
minimum if the non-minimal connection constant ^ > 0.123. Utilization of the Einstein equations at the minimum 
gives the stable configurations of the wormhole. For instance, in the case of conformal connection, ^ = 1/6, we found 
relation between the radius a of wormhole and mass m of the scalar field: am ~ 0.16. The Einstein equations say 
that the wormhole has a radius of throat a « 0.0141/pi and the mass of scalar field m « 11.35mpi. Therefore, this 
kind of wormhole, if it exists, may possess sub-Planckian radius of throat and it may be created by a massive scalar 
field with suppcr-Planckian mass. Obviously, the validity of the results obtained are restricted by the model taken - 
short-throat fiat-space wormhole. 

The goal of this paper is to consider the wormholes with more real geometry of throat and energy of quantum 
fluctuations of massive scalar field as a source of this background. The main problem in this case has rather math- 
ematical character. Even for simple profile of throat it becomes impossible to obtain a full set of solutions of radial 
equation in order to find the energy density of quantum fluctuations in close form. Nevertheless, it is possible to make 
some predictions about the existence of the wormholes by considering the heat kernel coefficients [il|. In fact, the 
crucial point is the existence of the negative minimum of the zero point energy. The sufficient condition for the zero 
point energy to have negative minimum is that the heat kernel coefficients B2 and B3 be positive |llj. This gives a 
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condition for parameters of model. More precisely, if a background is described by a parameter r with dimension of 
length and the domain where the space-time is "mainly" curved is defined by this parameter, then for small size of 
curved domain, r — s- 0, the zero point energy shows the following behavior 



327r2 

and in opposite limit r ^ c» we have 



If both these conditions are satisfied one can expect that the system will stay in minimum of energy which is 
characterized by specific values of parameters of wormhole and constant of non-minimal connection ^. Next step 
is the utilization of the Einstein's equations with energy-momentum tensor of quantum fluctuations as a source. 
Integration over volume the t — t component of this equation gives an additional relation between parameters of 
wormhole and zero point energy by using which we obtain the size of wormhole and mass of scalar field in terms of 
the Planck length and Planck mass correspondingly. At the beginning we may expect [TH that the size of wormhole 
and mass of field will be in the Planck scale. For this reason we are interested only in finding of the domain of the 
wormhole's parameters and constant non-minimal connection ^ for different models of the wormholc's profile. 

The manifest expression for coefficient B2 exists for arbitrary background, but this is not the case for coefficient 
B^. For this reason we adopt here the zeta-regularization approach (see Sec. IIIip . in frame of which it is possible 
to calculate the heat kernel coefficients and zero point energy itself. We pursue here another goal - to evolve the 
zeta-function approach for situations where it is impossible to find the full set of solutions of radial equation in closed 
form. We find a method to calculate the heat kernel coefficients in the background of wormhole with arbitrary profile 
of throat by using the WKB approach. Moreover, we obtain expressions for arbitrary heat kernel coefficients and we 
reproduce them in manifest form up to B3 for arbitrary profile of wormhole's throat. 

The organization of the paper is as follows. In Sec. ^we consider the geometry of wormhole with smooth throat. 
In Sec. mil we discuss the method of zeta-function for calculation of zero-point energy. The WKB approach for scalar 
massive field is considered in Sec lIVI The heat kernel coefficients are obtained in Sec0 We calculate them in manifest 
form for arbitrary profile of throat. The specific profiles of throat are investigated in Sec I VII and I VIII In Sec lVIIll we 
discuss the results obtained. Appendix A contains some technical formulas which are rather complicated to reproduce 
them in the text. 

We use units h = c = G = 1. The signature of the space-time, the sign of the Riemann and Ricci tensors, are the 
same as in the book by Hawking and Ellis 0- 

II. A TRAVERSABLE WORMHOLE WITH SMOOTH THROAT 

The metric of space-time of wormhole which is under consideration has the form below 

ds^ = -dt^ + dp^ + r^{p){d9^ + sin^ Odip^). (1) 

The radial variable p changes from —00 to -l-oo. In the paper we restrict ourself by wormholes with symmetric throat 
which means that r{—p) = r{+p). The radius a of throat is defined as follows a = r(0). We suppose that far from 
the wormhole's throat the space-time becomes Minkowskian, that is 

hm ^-^ . 1. 

p— >±oo p'^ 

The non-zero components of the Ricci tensor and scalar curvature have the following form 

2r" 



K = - 



n = - 



2(-l + r'2 + 2rr") 



The energy-momentum tensor corresponding to this metric has diagonal form from which we observe that the source 
of this metric possesses the following energy density and pressure: 

_ -l + r'2+2rr" 
^ " 8^ ' 
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-1 + r' 



In the paper we obtain general formulas for space-time ^ with arbitrary symmetric function r{p) obeying above 
Minkowskian condition. Two specific kinds of throat's profile will be considered. In the first model the profile of 
throat has the following form: 



r{p) = 77 + ^, (2) 

where a is radius of throat which characterizes wormhole's size. The embedding into the three dimensional Euclidean 
space of the section of the space-time by surface t = const, 6 = Tr/2 is plotted in Fig^I) for two different values of 
radius of throat. In Euclidean space with cylindrical coordinates (r, ip, z) this surface may be found in parametric 
form from relations: r = r(p), z'{p) = \/l — r'^ . In this background there is the only nonzero component of the Ricci 
tensor which reads 



" (p2 + a2)2- 

The second model has been considered in Ref.^^ and it is characterized by the following profile of throat 

r{p) = pcoth{-)-T + a. (3) 

T 

This model possesses a more reach structure. There are two parameters t and a. The latter parameter is radius 
of throat. In this model we may introduce another parameter which may be called length of the throat. The point 
is that the function r{p) turns into linear function of p starting from distance p > t/2 and the space-time becomes 
approximately Minkowskian. Therefore, the length of throat / = r. Using new variables y = p/a, a — r/a one 
rewrites the function r in the form below 

y 

r{y) — a{ycoth{ — ) ~ a + 1). 
a 

The parameter a is the ratio of the length and radius of throat. This parameter will play the main role in our analysis. 
It allows us to consider wormholes of different form, that is with different ratio of the radius and length of throat. 

In Fig^II-IV) the sections t — const, 6 ~ Tr/2 of this wormhole space-time are shown for different values of a and 
T. Namely in Fig^II) we represent two wormholes with the same radius of throat but with different length, and 
vice versa in Fig^II) we depict two wormholes with the same length of throat but with different radii of throat. In 
last picture Fig^IV) two wormholes with the same ratio of length and radius of throat, but with different values of 
throats' radii are depicted. Therefore the size of wormhole with the same ratio of length and radius throat is managed 
by parameter a. The parameter a describes wormhole's form. 



III. ZERO POINT ENERGY: ZETA-FUNCTION APPROACH 

We exploit the zeta function regularization approach 0, developed in Ref. and calculate the zero point energy 
of massive scalar field in this background with consequent using the Einstein equations. Let us repeat some main 
formulas from those papers. In framework of this approach the zero point energy 

^(^) = E E(^(n).. + ^^'y^'^' - r''Cc{s - 1), (4) 

3 (") 

of scalar massive field $ is expressed in terms of the zeta-function 

3 («) 

of the Laplace operator £ = —A + + SJZ. Here A = g'^^Vk^i is the three dimensional operator. The eigenvalues 
A(„) j + m? of operator C are found from boundary condition which looks as follows 

vI/(„)(A,i?) = 0, (6) 
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FIG. 1: First figure (I) represents the section t = const, 6 — n/2 of wormliole's space-time witli profile function r{p) = 
+ for two different values of radius of the throat. Three next figures illustrate the wormhole with profile of throat 
r{p) — /9coth(p/r) — t + a. Figure (II) and (III) illustrate that the a and r are radius and length of throat, accordingly. In 
last figure two wormholes with different a but with the same ratio radius and length of throat are depicted. It is seen that the 
parameter a characterizes the "size" of wormhole and a describes the "form" of wormhole. 
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where R denotes some boundary parameter. The solutions A = A(„) j of this equation depend on the numbers (n), and 
additionally they have the index j = 1, 2, . . . , which numerates the solutions of the boundary equation. Therefore, the 
zeta- function is a sum of expressions which depend on zeros of function ^'(n)- Next, according to Ref.j3| we convert 
the series over j in zeta function to integral and arrive to the formula 

2 n dk 

where the function ^ („) iii imaginary axes appears. 

The expression Q is divergent in the limit s — > we are interested in. For renormalization we subtract from E[s) 
all terms E'^™{s) which will survive in the limit m — > oo: 

£;*''(s) = lim E{s) 

m — ^oo 

and define the renormalized energy as follows: 

£:''™ = lim{E{s) ~ E'^'^'is)). (8) 

s—fO 

Because the pole structure of zeta-function does not depend on the value of parameters, it is obvious that in the limit 
TO ^ oo the divergent part will have the structure of DeWitt-Schwinger expansion, which has the following form 

X |so™*r(s - 2) + By2m^^is - |) + Bim^r{s - 1) + By^mTis - i) + B2r(s)| , 

where Ba are the heat kernel coefficients. In order to extract the divergent part of the energy we use the following 
procedure f3|. We subtract from and add to integrand the uniform expansion of In^E* up to to". We denote this 
expansion as (In ^E" Therefore, according to this, we represent the energy as the sum 

E{s)=Ef,n{s)+Eas{s) (10) 

of finite (in the limit s ^ 0) part 



1 /"OO 

EMs) = -i/^^^E^/ dk{k 

(n) 



2 - m2)i/2- 



X ^(ln*(„)(zA:,i?)-(ln*(„)(zA;,i?)r), (11) 
and the remains, which will be obtained from the uniform expansion part 

Basis) = -If^'' r " ^y^'-'i(^^'^(n)i^k, R)r. 

The last expression contains all terms which will survive in the limit m — > oo. 
Taking into account the obtained expressions in Eq. ((HJ we arrive at the formula 

E^'^"=Ef.,^ + E[r, (12a) 

where 



Ef^n = ^/<n(0) = -— E / dk^W^ 

, N Jin 



The divergent part E'^™ is given by Eq. @ 



(") 

X ^(ln*(„)(*fc,i?)-(lnvl/(„)(zfc,i?))'^^), (12b) 
= lim(i?a.(,s)-i?*''(s)). (12c) 
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The finite part Efi„ is calculated numerically. The second part, in practice, is found in the following way. By using 
the uniform expansion (In we calculate in manifest form the Eas{s) and after that we take the limit m — > oo in 

the expression obtained (the pole structure does not change) . All terms which will survive in this limit constitute the 
DeWitt-Schwinger expansion @ which we have to subtract in Eq. H12c|l . This way of calculation is more preferable 
because we may obtain the heat kernel coefficients in manifest form. The calculations of heat kernel coefficients in 
framework of this approach shows that the approach is suitable for both smooth background and for manifolds with 
singular surfaces codimension one fu] and two ^3 j the general formulas for which were obtained in Ref s. B and 0] ■ 

From consideration above we may find the zero-point energy for large and small size of wormhole Let the 

parameter a characterize the size of wormhole. In this case the E'^'^'^ /m is a dimensionless function and it depends 
on the parameter ma and some additional dimensionless parameters which characterize the form of wormhole. For 
example, in first model Q there is the only parameter a which is the radius of wormhole's throat and it characterizes 
at the same time the size of wormhole as a whole. Therefore in this model the E^^'^ /m depends on ma and there are 
no additional parameters. In the second model ^ there is an additional parameter a = t /a except parameter ma. 
For this reason the dependence of the zero point energy E'^^"/m on the mass is the same as parameter a. Because 
for renormalization we subtracted all terms of asymptotic over mass expansion up to B2 the asymptotic ma — > cx3 is 
the following 

m ~ 327r2m3 327r2(ma)3' ^ ' 

In opposite case, ma 0, the behavior of energy is defined by coefficient B2 (see also Eq. (0U): 

j^ren Inima)"^ ln(ma)2 ^ ^ 

m Azt: m ilir^ima) 

Here 63 and 62 are dimensionless heat kernel coefficients which may depend on the additional parameters. Therefore 
from these expressions we obtain the following sufficient condition that the zero point energy has a minimum: both 
B2 and B^ have to be positive. An additional condition may be obtained from Einstein's equations (see Sec lVIIVlTll . 

IV. MASSIVE SCALAR FIELD IN WORMHOLES BACKGROUND: THE WKB APPROACH 

We consider massive scalar quantum field in this backgrounds as a source for this space-time. In the frameworks of 
the approach used one has to find the spectrum of Laplace operator C: 

(-A + ^7^)$ = A$. 

Taking into account the spherical symmetry of problem we represent equation in the following form 

where YJ'"(0, (/?) are the spherical harmonics, / = 0,1,2,... and m e [—1,1]- The radial part of the wave function 
subjects for equation 

+ " - = -X% (15) 

To find the spectrum A we have to impose some appropriate boundary condition. It does not matter which kind of 
boundary condition will impose because in the end of calculation we will tend this boundary to infinity. We use the 
Dirichlet boundary condition in the spheres with radii R: p — ±R. For simplifying formulas we will work here with 
function ({s) = m'^'^Qcis)- With this notations the regularized ground state energy reads 



Because we need the solution for the imaginary energy only (see Eq. we change the integrand variable in 

radial equation H15|l to imaginary axis: A — > ivk and rescale for simplicity the radial variable: pk —> x. Therefore we 
arrive at the following equation [v — I + 1/2) 



+ 2-0-//2(l + ^)0+(-^ -e7^fc)0 = O, (16) 



7 



where the dot is the derivative with respect to x; = r'^k'^, and TZk = TZ/k"^. 

A general solution of radial equation (|16|) is the superposition of two linearly independent solutions 

'^{wk,x) = Ci(j)i{x) +C2M^)- 



(17) 



The first function (pi tends to infinity far from throat, for p oo, and the second one tends to zero. We consider the 
behavior of functions only for one part of space-time namely, with p > 0. The behavior of the solutions in the second 
part of space-time with negative p is found as continuation of the solutions from positive part of space-time. Now we 
impose the Dirichlet boundary condition at spheres p = ±R: 

W(zi/fc, +R) = CiM+R) + C2M+R) = 0, 
-R) = CiM-R) + C2M-R) = 0- 

The solution of this system exists if and only if the following condition is satisfied 

^liii^k, R) = M+R)M-R) - M-R)M+R) = O- (i8) 

The contribution from the second term in equation above is exponentially small comparing with first one in the 
limit R ^ 00. In order to see this let us find the uniform expansion of solutions 0i and (/)2. Moreover, we need this 
expansion for renormalization and calculation of the heat kernel coefficients. Let us represent a solution in the 
exponential form 



1 



'lav 



(19) 



where a = r(0), and substitute it in the radial equation (|16() . One obtains a non- linear equation 

0. 



We represent now the solution in the WKB expansion form 



00 



and substitute it in equation above. This gives the following chain of equations 

^-1 = 



5o 
^1 

Sn+l 




SkSn — k + 2 — Sn 

rk 



1,2,- 



(20a) 

(20b) 
(20c) 

(20d) 



There are two solutions of this chain corresponding to sign in the first equation. The sign plus gives the growing 
(for positive coordinate p) solution which we mark " -t-" and sign minus gives solution which tends to zero at infinity 
which we mark by sign "— ". Therefore 



M+R)M-R) 



,S+{+R)+S-(-R) 



lav 



To find an expansion for the sum ^^(-l-i?) + S (— -R) we need the following properties of function >S'^(a;) 

^InK^) = +'5'2n(2;), 
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and 

where n = 0, 1, . . . . The first two equations are the consequence of the structure of chain and the last two equations 
are due to the symmetry of metric function rk{x) = 7fc(— x). 
Taking into account these properties of symmetry we have 



S+{+x) + S-{-x) - Y."^'^" 

oc 



r!=0 



(21) 



S+{+x) + S-{+x) ^ + c:,^-i] 



(22) 



-2n 



+ C'jn +'2- I Stndx 



+x 



Here the C„ are the constant of integration of system 12U|I . 

Therefore we may express the combination we need (|21|) in term of the H22() : 



S+{+x) + S-{-x) = S+{+x) + S-{+x) + Y r S^n-l 

n=0 "'-^ 



dx 



To find the combination S'+(+a;) + S {+x) we use the Wronskian condition. Because these solutions are independent, 
they obey to equation {ak = ak) 



w{M^),Mx)) = -. 



The origin of this relation is the following. Suppose we try to find the scalar Green function of the Klein-Gordon 
equation: 

5^{x,x') 



(S^'^V^V, - - e7^) G(x, x') = 

V-5(a;) 

in background (Q) . It is very easy to extract time and angular dependence of the Green function 



(23) 



/ + 00 7 ^ ^ 

?E E yr{e,v)Y-^{e',v')e 



;=0 ra=~l 



and we arrive to equation for radial part of Green function which reads (A^ — llP' 



-e7e|(/.(p,p') = 



Kp-p') 



or in dimensionless variables (A — > ivk, kp — > x) 



, k5{x — x) 



di + 2^5, - v\\ + ) + (-^ - e7^fe) 0(x, x') 

As usual we represent the radial Green function in standard form: 

<j){x,x') = 9{x' - x)<j)i{x)(t)2ix') + 9{x - x')02(a;)0i(x'), 

where (j^i and 02 are two linearly independent solutions of homogenous equation and (/>i tends to infinity for p — s- oo 
and 02 tends to 0. The Wronskian condition appears if we substitute this form of radial Green function to radial 
equation above: 



VF(0i(x),02(x')) 
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Therefore, if two functions 4>i and 02 describe the system, they have to obey this Wronskian condition. 
For solution in exponential form (|19l) this condition gives 

gS+(x)+S-(2:) ^ ^i^flfc 1 

S+ix)-S-{xy 

The denominator in rhs has the fonowing form 

CXD 

n=0 

Taking into account these two expressions above we arrive at the formula 

S+{x) + S-i^x) = Ha,) - HSUr,^) + J] / S+^_,dx - In 1 + ^ ^"'"^ • (24) 

^ n=0 "'-^ I n=l -^-l J 

The main achievements and peculiarities of this expression are: (i) the rhs is expressed in terms of derivative of 
functions 5+, we do not need to find the constants of integration in the chain of equations if^ . (ii) the odd and even 
power of v are separated, which leads to separation of the contribution to heat kernel coefficients with integer and 
half-integer indices, (iii) the rhs is expressed in terms of functions Sn with odd indices, only. The first three functions 
S2n-i are listed in Appendix^ formula ljAla|) . We would like to note that this formula is valid for arbitrary, but 
symmetric, r(p) = r(— p), metric coefficient. 

From this expression we may conclude that the contribution from the second term in condition (|18() is exponentially 
small comparing with first one. Indeed, the main WKB term in Eq. (|24|l gives the following contribution: 

k 1 f /■+''^ ■ 

2^ S-ir, [ J-kR 

k 1 f f+^^ 1 

0i(-i?)02(+i?)«:^^^exp<^-i. / St^dx). 

^"SZir^' [ J-kR J 

Because the function S^^ is positive for arbitrary R we observe that the second expression gives exponentially small 
(for R — !■ oo) contribution comparing with first one and we will omit it in what follows. 



V. HEAT KERNEL COEFFICIENTS 

Let us now proceed to evaluation of the heat kernel coefficients (HKC). The formula (|24|l allows us to find HKC in 
general form for arbitrary indices. Taking into account above discussions we have the following expression for zeta- 
function 

C(. - i) = Y: / dkik' - ^ f"-'^^ {S^{+R) + S-{-R)} . (25) 

To find the heat kernel coefficients we use the uniform expansion given by Eq. 124|) . As it will be clear later, the 
odd powers of v will give contribution to HKC with integer indices and even powers of v produce the contribution 
to HKC with half-integer indices. The well-known asymptotic expansion of zeta-function in three dimensions has the 
form below 

Casis - i) = (J)3/2r(,-i) E {-'-"5^r(. + / - 2) + m--s,^,/2r(. + ^ - ^)} • (26) 

For simplicity we introduce the density of HKC with integer indices J5; by relation 

n+R _ 

Bi = / dpBi{p) 

J-R 

and first of all we will obtain formulas for this density. 
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Let us consider the part of Eq. (|24(l with odd degree of v. The contribution to zeta-function is the following: 



1, o,. 2cos7rs 



We change now the variable of integration x — kp and take the derivative with respect to k: 



kR 



cr(--i) = -m2^±^E-'"'^ / dkkik'-'^f/'-^Y.-"'"' / -2p-i(fc,p)dp. 



The first four functions S2p-i are listed in Appendix 1X1 formula l|A2a|l . The general structure of these functions is 
the following 

2p 

S2p-i = a2p-i,„^:~^~"~^, 

where a2p-i,n are the functions of r{p), and z = 1 + k'^r^{p). 
Next, we integrate over k using the formula 

dfcfc(fc2 _ + er^yi = 1^2.-3^-3+2.+2,(^2 ^ ^2^2) j -.-g r( § ~ ■s)r(g - | + ,9) 

i/^ 2 r(g) 

and obtain the following expression for odd part of the zeta function 

1 „,2s f+B rfq + Ti + n-n ,,2«+l 

2 r{s-^)J_R f^f^o^O r(j3 + n+i) (jy2+„^2^2)6+p+« 1 

By using the binomial of Newton we reduce the power of v in nominator 

V = - Vf-mVr-'^^^ Z{s+p + n-l-q) 

Z-.^ (y-^ + 'm?r-^Y+P+^-^ 2^^ ' q\{n-q)\T{s+p + n-l-q)' ^ ' 

where 



z(.) = 2r(.)E 



1=0 ^ ' 

To obtain the HKC we need asymptotic (over mass m) expansion of the zeta-function. The asymptotic expansion 
of function Z(s) was obtained in the Ref. and it has the form below 

oo 

Z{s) = (TOr)-2^ E ^i(s)(mr)-2', (29a) 

= r(.s-l), (29b) 
Ai{s) = 2^^r(; + ,s)Cff(-l-2?,i), (29c) 

where the C^ii{a,\)) is the Hurwitz zeta-function. 

Taking into account formulas (|28(l and H29(l one has the following asymptotic series for odd part of zeta-function 

1 1 /• + -R oo ; 2p n -r.lt i i \ 

crc-i) ^ ^.i: £ E i: +; r;' 

n\ „_^ A;_p- i(s+p + n- 1 -g) 



(7!(n — gr)! r(s + p + n - 1 - g) 
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As it was expected at the beginning this is series over even degree of mass and it gives contribution to the HKC 
with integer indices. Comparing above expression with general l|26|l we obtain general formula for arbitrary HKC 
coefficient with integer index 

Therefore, to obtain HKC with index I we have to take into account expansion up to z/^"^'. 

Let us now proceed to HKC with half-integer indices. To find them we have to take into account the rest part of 
Eq. ^ with even powers of v in expression for zeta- function (|25|l . The general form of even part is 

1 ( °° o+ ] °° 

(S+ix) + = In(afe) - - HS\r,f) - In 1 + J] ^-^"^ = ^ ^-^pe,,, (31) 

L ri=l —1 ) p=0 

where the first four functions i?2p are listed in Appendix IXI formula IjAlbjl . 
We substitute now the expansion (|31|l in the expression for zeta-function: 



and take the derivative with respect to the k: 

Cris -\)= ^m^^t^ I dkk{k' - ^y/'-^ E .~'^s,,. (32) 

i=0 •'17 ^ p=0 

The functions S2p have the following structure 

2p 

■S2p = E/ ^'^V,riZ P n 1^ 

where z — 1 + k'^r^{R). The first three coefficients are listed in manifest form in Appendix A (see Eq. (|A2b|) '). The 
coefficients a2p,„ depend on the parameters R and a and do not depend on the variable of integration k. Going the 
same way as we did for HKC with integer indices we obtain the following asymptotic expression for even part of the 
zeta-function: 

^''^ (^-2^ " m73T)Z.Z.Z.Z."2p.nm r ^(p + n+l) 



n 



g!(n-g)!^ T{s+p + n-^-q) 

As it was expected the even part of zeta-function is the series over odd powers of mass and, therefore, it gives 
contributions to HKC with half-integer indices. Comparing this expression with general asymptotic series for zeta- 
function we obtain the following formula for HKC with half- integer indices: 

4^3/2 ^,, r(p + n - 1 + .s) n!(-l)"-^ Ai_.p_,{s+p + n - ^ - q) 

'^'^~ns + l-l)}-,hh"""-''' r(p + n + l) ,!(n-,)! T{s+p + n-\-q) ' ^''^ 

We would like to note that the right hand side of formulas (|30|l and H33(l does not depend, in fact, on the s, which is 
confirmed by straightforward calculations. 

These formulas look very complicate but calculation may be done easily using simple program in package Mathe- 
matica. Indeed, the functions ^^(a;) and i?2n may be found by using formulas H2()|l and H31|l. The functions s„ are 
obtained from the following relations: 



= kdk 



x—kp 



1 d 



12 



The first four HKC coefficient (density) with integer indices are Usted below 



Bn = 47rr^ 



Bi = 

= 



»7r 

Y 

Sire •■ 

315r2 
3r4 



r'' + rr" 



8^ U - 



1 



-1 + r' +2rr" 



6 

2 87r£ 
3r^ 



-l + r' +2rr" 



(-1 + r'2) - r (-5 + 7r'2) r" + Tr^rV'^) + SrV^^) 
2 (-21 + 17r'^) - 6r (-35 + 59r'^) r" + 21r2 (^Tr"^ + 2Ar'r^^^^ + 210rV('') 



(34a) 
(34b) 

(34c) 



-l + r' +2rr" 



3r4 



(-1 + 



(1 + 9r'2) - 2r (-1 + r'^) (-5 + Or'^) r" 



47r^ 
315r4 



[-42 (-1 + r'2) (1 + 15r'2) + 2r (-252 + 105/^ + 859/^) ' 



- 2r2 (^-525r"2 + 2517/^"^ - 420rV(^) + 808/^(3)) + ir^ {iO^r"^ + 1354rV"r(3) - llbA^^ + 271r'V('')) 
+ 21r4 (27r(3)2 + 27r"r(4) - 13rV(5)) - lOSrV^^) 



45045r 



J [4 (-572 - 9009r'2 + 9341r" 



4r (-6006 - 15015r'2 + 62039r"*) r" + 13r2 (^-4620/'^ + 32943/^"^ - 4620rV(3) + 11564r'V(^)) (34d) 
286r3 i 308r"^ + 1139rV'V(3) - lOSr^") + 22?,r'^r^^A - 429r4 ( 47r(3)2 + 24r"r('*) - lAr'r'^^A + 12012rV(s) 



In above formulas the function r depends on the radial coordinate p whereas the heat kernel coefficients with half- 
integer indices 



Bl,2 = -4^3/2^2^ 



B 



3/2 



3/2 



B5/2 - 



+ 



-871^/2^ [-1 + r'2 + 2rr"] + ^ [-4 + 3r'2 + 6rr"] 



(35a) 
(35b) 



-8^3/2^2 



7r3/2 
120r2 



[-l + r'2 + 2rr"]'- 



2^3/2^ 



3r2 



4 (-1 + r'2) - lOr (-2 + 3r'2) r" - 3r2 (4r"2 - 3r'r(3)) + GrV^^^ 



2 (-16 + 15r'4) - 5r (-32 + 63r'2) r" - lOr^ {br""^ - 14rV(3)) + 90rV('^) 



(35c) 



depend on the radial function r at boundary: r = r{R). From Eqs. (|34|l and (|35|l we observe that the HKC Bi and 
Bi+i/2 are polynomial in ^ with degree I. 

It is well-known 0| that the heat kernel coefficients with integer indices consist of two parts. First part is an 
integral over volume and another one is an integral over boundary. We obtained slightly different representation for 
this coefficient as integral over p. But it is easy to see that they are in agreement. Indeed, let us consider for example 
coefficient Bi. According with standard formula wc have 



TZdV + - I trKdS 
V 3 



p=-R' 



Volume contribution is exactly the same as we have already obtained (|34b|) . Surface contribution from above formula 
is 



3^ 



trKdS 



1 



IS 3 
From our result Ij34b|l we get the same expression 



- / trKdS 



p=-R 



87r 

T 



-R 



+R 



167r , 

r r 

3 



167r 



p=R' 



dp = — / [rr']' dp = —z—rr' 



-R 



p=R' 



It is not so difficult to verify that the heat kernel coefficients up to B2 are in agreement with general expressions. 
There is no general expressions for higher coefficients. 
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According with Ref. 'll'| the sufficient condition for existing of the the self-consistent wormholes may be formulated 
in terms of two heat kernel coefficients 

/ + 00 
-oo 

/ + 00 
dpBs ^ 713,3^^ + h3^2C^ + hs^iC + h3,0- 
-oo 

Namely, both B2 and have to be positive The coefficients hk^i of polynomials depend on the structure of 

wormhole. Therefore the problem reduces to analysis of polynomial in of second and third degree, the coefficients 
of which depend on the structure of wormhole's space-time. Wormholes with different forms may exist for different 
values of non- minimal connection ^, or vice versa for some ^ the above polynomials will be positive for specific form 
of wormholes. 



VI. THE MODEL OF THROAT: r(p) = ^ + 

In this section we consider in detail the specific model of wormhole with the following profile of throat r(/9) = 
a/p^ + o?. From general expressions (|34l) we obtain the density of heat kernel coefficients with integer indices which 
are 

Bo = 47rr^ 

B3 = (2583561a6 - 31574380^^ + 751820aV - 480rS) -|- ^„ (47263a* - hl^Ma^r^ 13540r*) ^ 

45045ri" ^ ' 315ri" ^ t y s 

- 3;acr(73«'-62r2)C'-H^^e' 

Integrating over p from — i? to +i? we obtain HKC. Below we reproduce their expansions in the limit i? ^ 00 up to 
terms 1/i?: 

Bo = 5^+87ra2i?, (36a) 
„ 167ri? „ 9 Is 167ra^ ,^ 1, n 

TT^ 327r 

S2 « - m + 3) - — (36c) 

20a 315it 

B3 w ^(5880^^- 6300^2^ 2226^- 257). (36d) 

4032a-^ 

The formulas for first three coefficients with half integer indices may be found from general expression H35() . Below 
we have listed them with their expansions for large value of R 



Bi = 


_4^3/2^2 


2 




B3 « 


^3/2 




2 


3 


Bs « 


^3/2 




2 


60i?2' 



-47r3/2(i?2 



60i?2 



Let us now proceed to renormalization and calculation of the zero point energy. As noted in Sec. IIIII fsee Eq. ffT^ ) 
we have to subtract all terms which will survive in the limit m 00. According with general asymptotic structure of 
the zeta- function given by Eq. 1)26(1 in this limit the first five terms survive, namely HKC up to i?2- Because the zero 
point energy is proportional to zeta-function we may speak about renormalization of the zeta function. According 
with Eq. (|12|l we take asymptotic expansion for zeta-function up to (in the limit m — > 00 these terms give 
asymptotic (over m) expansion (|26|l up to heat kernel coefficient B2) and subtract its expansion over m up to B2 
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from it. After taking the limit s ^ we observe that this difference will give -E^™ (|12c|l . First of all we consider this 
part and later will simplify the finite part (|12bp . 

We should like to make a comment. In the problem under consideration we have two different scales: R and a which 
give us two dimensionless parameters mR and ma. To extract terms for renormalization we turn mass to infinity, 
which means the Compone wavelength of scalar boson turns to zero and becomes smaller then all scales of model. In 
other words it means that we turn to infinity both parameters mR and ma. After renormalization we will turn mR 
to infinity separately in order to obtain the part which does not depend on the boundary. 

Let us consider separately two parts of asymptotic expansion of zeta-function according with odd and even powers 
of h'. First of all we consider odd part which gives the HKC with integer indices. All singularities are contained in the 
first three terms in Ed. H2t)|l with Bq, Bi, i?2- After subtracting these singularities we tend s — > and obtain some 
infinite power series over parameters mp and ma. Next, we have to integrate over p and tend mR —^ oo. For this 
reason we have to obtain some expression instead of series to take this limit easily. It is impossible to take this limit 
directly in series. We will use the Abel-Plana formula to extract the main contribution from series in this limit. The 
rest will be a good expression for numerical calculations. Moreover, from this rest part we will extract terms which 
will be divergent in the limit ma to analyze our formulas. 

Our starting formula for odd part of zeta-function is Eq.(|23) which cut up to p = 2: 

Expanding the denominator with the help of the formula 

n I ,2.-. f^ (-i)"r(n + .) 

we represent Eq. H37|) in the following form 

1 1 1 r + R °° 2 

& - -) ^ (,^^3/.r(.-l) E E-'"/n,p(^), (38) 



where 



In order to make formulas more readable we make everything dimensionless but save the same notations. In any 
moment we may repair dimensional parameters by changing R mR an a ^ ma. In this case we rewrite the 
expression for zeta-function (|37|l in the following form 



1 m m^" /"^^ °° ^ 

^-'^^-2)=(4.)3/2r(.-i)L ^^EE^-(^)' 



(39) 



From this expression we observe that for p — Q the first three terms are divergent with n = 0, 1, 2; for p — I the 
first two terms with n = 0, 1 and at last for p = 2 the only term is divergent with n = 0. We remind that 

r(s - n)_^ = tip. Q + ^,(^n + 1)^ + Ois), CHis + 1, g),^o = ^ - *(g) + 0(s). 

For this reason we represent zeta-function H37() in the following form (for s — > 0) 

Ct^(.s-i) - (4 "3/2 rr. '^^]-"E/-o(.s) + E/'.o(0)j (40a) 

r^^E/'M(.'^) + E/".i(0) (40b) 

n=2 J 

/n,2(0) I (40c) 



(47r)3/2 r(s - 




m 1 




(47r)3/2 r(s - 




m 1 





(47r)3/2 r(s - i) 
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and will analyze each part separately. 

To illustrate the calculations we consider in details the first part 14Ua|) . First of all it is not difficult to find the 
manifest form of singular part in the limit s 0: 



n=0 



47rr2r(s - 2) + -r(s - 1) 
o 



120r2 



r(s) +7rr2[-3 + 47 + 81n(2)] 



-[1 + 21n(2) + 24C^,(-1)] - T^[-7 + 21n(2) - 1680C^^(-3)]. 



120r 

We observe that this term gives contribution to Bq-, Bi, and B2 according with gamma functions. For renormalization 
we have to subtract from this expression the first three terms according to our scheme. 

There is one important moment which is crucial for our analysis. Above formula contains all terms which survive 
in the limit s — *■ for arbitrary mass of field. For renormalization we have to subtract asymptotic expansion in the 
form (|26l) . There is a difference in factor r^". For this reason after renormalization factor 



(r^^ - l)(47rr^r(s - 2) + -T{s - 1) 



Ttt 



■r(-)). 



27rr" — — 



Ttt 
120r2 



Inr^ 



3^ ^' ' 120r2 

appears. If we take into account all terms in Eq. (|40|l we obtain the following contribution 

ln(r2)(lso-Bi+B2). 

Exactly the same structure was observed before in Refs. ^| . This term defines the behavior of energy for small 
size of wormhole because it is maximally divergent for small size of wormhole. 
Therefore the renormalized contribution is 



(41) 



r!=0 



27rr^ 



Ttt 
120r2 



In + nr^ [-3 + 47 + 8 ln(2)] 



-[l + 21n(2) + 24C^(-l)]-^ 



[-7 + 21n(2) - 1680CflX-3)]- 



We represent the finite part in the following form 



n=3 1=0 ^ ^ 



ill 

2^ 



In 1 



by using standard series representation for Hurwitz zeta-function. To find more suitable form for these series we use 
the Abel-Plana formula and obtain 



E-^ 



1=0 



In 1 





to 






1 




2 





-r^ ln(r2) - r2[21n(2) + j - -] + 2 T In 



2 



lr4ln(r2)-lr2+r4[21n(2)+7. 



dvv^ 



■In 



Taking into account these formulas we have 

2 00 
E/«^o(0) + E/ro"(0) = 167r 



n— 3 



dvv'' 



1 



In 



9 

1 



^2 1^ 



We now integrate this formula over p from —R to +_R according with Eq. (|40al) and take the limit R 
this we arrive at the expression 



m 
~167r2 
77r2 



/a 



m 

"l67r2 
a / 7 



where 



a 



1 



^ 60^ + - ( + T7^ 1^(2) + 14C^(-3) 



120 10 



After 



The manifest form of the regular contribution has written out in Appendix (see Eq. IA3(I . We extracted all terms 
with logarithm and that which is singular for a — and collected them in Z^™^. The rest part, Wa, is a regular 
contribution. 
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FIG. 2: The plot of the summary contributions: f ~ fa + ft + fc and ljs = i-^a + i-^b + u>c for £, ~ \: a) summary contribution 
/ and b) regular part lus 



Using the same procedure for second and third parts we obtain the following expressions 



167r2 



rsinq 

fb 



IGtt 
1 

6 

m 

167r' 



1 /2 



ln(a) + - 



-l(l + 24C^(-l))e+i(l + 18C^,(-l)) 



6r - 



20 



where 

ln(a) 



TT 

a 



i(-7 + 127 + 361n2)e^ + i(15 - I67 - 481n2)^ + T^i-U^ + I267 + 378 In 2) 
2 D 6W 



The manifest form of the regular contributions lui, and cUc have written out in Appendix 1X1 fEo. (jA3|) ). 
Putting together all contributions in (|4()|l we obtain 



Co. 



ren 
dd 



1677^ 



\n{o?y ( - 
' a 



where 



to = + CJ6 + CJc + 



i(-7 + 127 + 361n2)C^ - i(-13 + 48Cfl(-l) + I67 + 481n2) 



180 



-41 + 360Ci7,(-l) + 2520C^j(-3) + 637 + 207 ln2) 



In Fig. 121 we reproduce a plot of sum of all three contributions: / = /„ + ft + fc, = + '^b + for ^ \- 

Let us now proceed to the contribution from even part of zeta-function. We start from Eq. H32|) and do not take 
the limit of great mass. Integrating over k we obtain the following expression for this even part 



1 2s °° 1 2p 

^even/ _ j_^ _ 

O' T( „ _ l\ 



a2p,k'> 



„2s-3 



r(s + p + /c- i) 



2k+l 



_g T{p + k+l) {,y2 j^^2^2y+p+k-y 

Here r is taken at the point R. Now we pass to dimensionless variables (with the same notations) and have 



(42) 



1 TD 

/■even/ _ _\ _ 

°" ^ 2' f47r)3/2rfs - i) ^ ^ 



1 ^ (47r)3/2r(,s+p + fc-l) 



r(p + fc + i) 



(43) 



where 



Pp.k 



E 



,2fe+l 



^ (,^2 ^.^2j,2)s+p+fe-i 
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Analytical continuation s ^ in this series may be easily done by using the Abel-Plana formula: 

1,1 f°° , , 2 r v'^dv 

00,0 = Tr:r -2r / dwEx{v) - - / , Ex{v), 

3 24 J, rj^ ^^^^ 

2 r I 

^1^0 = -r+-J diyiy[3Ex{iy) +vEx'{i^)]^l- ^ 



91,1 



-2r+- / dvv\&Ex(v)+QvEx'{v) + v^Ex"{v)]\ l- — 
r Jo V r-^ 

'i 2 r 1 I 

--r+- / diyvlSExM + \2vEx'(v) + Av'^Ex"{v) + -v^Ex"'{v)]\ 1 
3 r Jo 3 V 



where 



Ex{v) 



Now we substitute these formulas in Eq. (|43|l and take the limit i? — s- oo. In this limit all integrals in expressions 
for 0p jj are smaller the l/r. Taking into account that ao^o ~ i" and ai.fc ~ we obtain in this limit (we repair 
dimensional variables) 

Therefore after renormalization (subtracting these two terms) we take the limit R oo and obtain zero contribution 
from this even part 

Thus there is the only contribution from the odd part. Collecting all terms together we arrive at the following 
expression for zero point energy 



pren _ ) ^2_2 ' ^ 



3C - £ + — 



8ma[^-^]]+n}, (44) 



where 



n = L, + i2nY,v^ I dy{y^--f'^ (45) 
1=0 ''i ^ 

X ^ (5+(+mi?) + S-{-mR) - \S+{+mR) + 5-(-7ni?)]"™""'''| 

and fc = my. 

The main problem now is the calculation of the last term in expression for fJ. Let us simplify the expression and 
show that in the limit i? — > oo the divergent parts are cancelled. Indeed, let us consider the first five terms of uniform 
expansion 



2 i.+kR 1 

[S+{R) + ^-(-i?)]™'""^- = / ^2+.-W2; + E 

n=0 J -kR 



It is very easy to take the limit in the part with even power of v by using the manifest form of E2n listed in Appendix 
IXI The only term which gives the non-zero contribution is Eq: 

-EoU^oo = ~2\n{kR) + \n{ka) + 0{R-^), 

E2\r-.oo = 0(i?-2). 

The part with odd power of u in the uniform expansion brings the single linear on R divergent contribution coming 
from S^i'. 

2{kR)v. 
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Therefore in the Umit i? — > oo the uniform expansion gives the foUowing divergent contribution: 

2{kR)iy - 21n(fci?) + In(fca). 

Because later we have to take the derivative with respect to k we may rewrite this expression in the foUowing way 

2{kR)iy-\n{ka). 

To take the hmit of large box in S^{+R) + {—R) let us reduce the radial equation to standard form of scattering 
problem by changing the form of radial function 4> rp/r{p). In this case the equation reads 



^p = X^ip. (46) 



This equation looks similar to the equation of scattering problem in one dimension (do not forget that p € 
{—oo, +00)) with non-singular symmetric potential 

/(^ + 1) , a^(l~2g) 
+ (p^ + a-^)^ 

From standard theory of one-dimensional scattering we know that there are two independent solutions which have 
the following properties (as opposite to traditional notations (t>{x) we use solutions (j){—x) to make coincidence with 
functions we use) 

^^^P> ~ \ e-*^'' + si2(A)e'^'', p->-oo' 

... r S2i(A)e-*-^'' + e*^'', p~*+oo 
"^^^P' I S22(A)e»^'', 

where Saa{X) constitute the s matrix of scattering problem. Due to symmetry of the potential the components of the 
matrix obey the relation S22 = sii- 

Now we change energy to imaginary axis: A ikv and obtain 

M+R)U-R)n^^ = s\Sku)e''""'- 

Therefore 

[S+{^R) + S-{-R)\ ^^^=\n{<l),{+R)<l)2{-R))j,^^=HsUikv)]+2kRv, 
and the divergent parts in (|46|) are cancelled. Thus we arrive at the following expression for 

1 ^l/2 



^ />00 1 

= L0 + i2nY,v^ I dy{y^-—Y 
1=0 •'i ^ 



X — <jln[ys?i(*yj.m)]-^^^j.i-2"y \St^_^- E^.^\dx\ (48) 



Thus, we express the finite part of the zero point energy in terms of the s matrix of scattering problem, namely in 
term of the transmission coefficient of the barrier in imaginary axis. Similar relation was found by Bordag in Ref.j2|- 
The potential V/(p, a) of scattering problem has the following properties 

V/(0, a) = ^ ^ — the height of barrier, 

/+00 7t\21(1 + 1) + 1 2£1 
Vi[p,a)dp — the work against potential barrier. (49) 
-00 2a 

Therefore the zero point energy has the form (|44|l . where the function is given by expression 148() . We should 
like to note that according with jj^] the factor before logarithm term in (|44|l is {B2 — i3i)fl^oo- The origin of this 
structure has been already noted in Eq. H41I) . 
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Now we analyze qualitatively without exact numerical calculations the behavior of energy for small and large radii 
of throat. According with Eqs. 1)14(1 and (|13|) the zero point energy in 3 + 1 dimensions has the following behavior 
for small and large values of throat 

^ _B2Harnf^ ^^^^ 



327r2 



or in manifest form: 



ln(r 



32 
m 



ma 
1 



327r2TO2 



, a — > oo, 



32 (ma)3 4032 



[5880^^ - 6300^2 + 2226$ - 257] , a ^ oo. 



(50) 
(51) 



It is easy to verify that coefficient after logarithm in Eq. (|5()ll . which is contribution from B2 in the limit R ^ oo, is 
never to be zero or negative. It is always positive. For this reason the zero point energy is positive for small radius 
of throat for arbitrary constant of non-conformal coupling $. In the domain of large radius of throat the expression 
in brackets in Eq. (|51|l . which is the contribution from in the limit i? — > 00, may change its sign. It is positive for 
f > 0.266 (energy negative) and negative (energy positive) in the opposite case. Therefore we conclude that there is 
minimum of ground state energy if constant ^ > 0.266. The situation is opposite to that which appeared in our last 
paper llj, where the energy for large value of throat (which was defined by -B5/2) was always positive, but for small 
radius of throat it could change its sign. 

Let us now consider the semiclassical Einstein equations: 



87rG 



til// 



(52) 



where G^^ is the Einstein tensor, and (T^,y) is the renormalized vacuum expectation values of the stress-energy tensor 
of the scalar field. The total energy in a static space-time is given by 



E = [ ey/^d^x, 
Jv 



where e 



-in 



spherically symmetric metric we obtain 

E = - 

The zero point energy has the following form 



Glc /SnG is energy density, and the integral is calculated over the whole space. In the 



,4 /"OO 



c 
2G 



Glr'{p)dp 



4 

c TTa 

^IG' 



he 



-f{ma,^). 



(53) 



(54) 



In the self-consistent case the total energy must coincide with the ground state energy of the scalar field. Equating 
Eqs. l(13|) and gives 

C^Qi Tic 

— = -/(ma,0, 



a ^ lp\/2f{ma,S,). 

Considering this equation at the minimum of zero point energy we obtain some value of wormhole's radius. The 
concrete value of radius may be found from exact numerical calculation of the zero point energy as function of ma. 
But without this calculation we conclude that the wormholes with the throat's profile ^ may exist for ^ > 0.266. 
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FIG. 3: The discriminant A — 62.1 ^ 462, 262,0 of polynomial in ^ as function of a. It is always negative for all values of a. It 
means that the zero point energy is always positive for small value of radius of throat. 

VII. THE MODEL OF THROAT: r{p) = pcoth f - r + a 

We will not reproduce here the density for heat kernel coefRcients in manifest form due to their complexity. They 
may be found from general formulas H34() . There are two parameters in this model a and a = r/a. The dimensional 
parameter a characterizes this kind of wormhole as a whole. Small value of this parameter indicates small size of 
wormhole. The dimensionless parameter a characterizes the form of wormhole - its ratio of the length and radius 
of throat. By changing the integration variable p = xa we observe that coefficients B2 and B3 have the following 
structure which is clear from dimensional consideration: 



— 1 

B2 ^ dpB2 = - [62,2^^ + 62,1c + 62,0] , 

J —00 ^ 

B3 = dpBs = — [63,3^=^ + 63,2^^ + 63,1c + 63,0] 



where h^ i depend on the a — r/a, only. We note that 62.2 > as it is seen from (|34c|l . Therefore we may analyze 
the zero point energy for different values of the parameter a. From the general point of view we have the following 
behavior of the zero point energy for small size of wormhole that is for small value of a — > 0: 

In(ma)^ In(ma)^ ^ 2 c^h 1 

Using the general expression for coefficient B2 it is possible to find in manifest form the polynomial in ^ in H55|l for 
great value of a — > cxd (small radius of wormhole throat comparing with its length): 

E-^-^^^^^^m'-m + i], (56) 

and for small value of a — > (small length of wormhole throat comparing with its radius): 

i?-" « _(15-^^)ln(ma)^ . 2 _ gQ^ + 7] . (57) 
13507raa L s j y j 

The numerical calculations of the discriminant A = 62 1 ~ 462,262,0 of polynomial in ^ as function of a is shown in Fig. 
13 From this figure and Eqs. H56() . H57|l we conclude that the discriminant is always negative for arbitrary values of a. 
It means that the zero point energy is always positive for small size of wormhole for arbitrary constant of non-minimal 
coupling C and arbitrary ratio of the length of throat and its radius. 

The behavior of zero point energy for large size of wormhole (a — s- cx)) has the following form 

« -^JT^ = - 00 2 2 3 [^3,3C' + 63,2c' + 63,1c + 63,0] . (58) 
327r^m^ 327r^m^a'^ 

The zero point energy will get minimum for some value of a if above expression will be negative. Let us consider the 
polynomial 

P = 63,3c' + 63,2c' + 63,1c + 63,0 (59) 
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FIG. 4: The plots of the polynomial Pa^ for different values of ratio a — t /a. 



for different values of a starting from small value of it. Zero point energy will get minimum if this polynomials is 
positive. In the limit a ^ we have approximately 

3^ 5127r (45 ~ 47r2) / 2567r (21 - 27r2) ^ 2567r (45 - 47r2) V 2 



135 " \ 189 

' 5127r (21 - 27r2) 2247r (45 - A-k"^] 



945 



675 



135 

\ 3687r(21-27r2) 447r (45 - 47r2) 
■a U + -a. 



6615 



2025 



In this expression we saved terms up to a 
real: 



This polynomial in the limit a ^ has two complex roots and one is 



5(21-27r2) 1 



^ 14 (45 - 47r2) a' 
Because the coefficient with is positive the polynomials will be positive for all 



5(21-27r2) 1 
14 (45 - 47r2) a' 



Therefore for small values of a = r/a we have low boundary for parameter ^ where the wormhole may exist (see 
FigEII)). The greater a the smaller the low boundary of For a > 1.136 the conformal connection ^ — 1/6 will be 
greater then the low boundary. At the point a — 1.26 two domains appear where polynomial is positive. First domain 
is 0.188 < ^ < 0.841 and second £_ > 0.841 (see Fig^lll) for a > 1.26). The low boundary of second domain will 
increase for greater a and it disappears for a = 1.65. At this point the coefficients 63^3 = and the polynomial turns 
out to be parabola (see Fig^III)) with positive part in domain: —0.088 < ^ < 0.358. For greater a we obtain upper 
boundary of ^ where the polynomials is positive because the coefficient with is negative. Starting from a = 1.65 
we have two domains where the polynomial is positive (see FigEJlV)). First one closes to —0.088 < £, < 0.358 and 
another one is smaller then some negative value of ^. For a — 2.08 the high boundary of second domain will coincide 
with low boundary of first domain and we get the only domain where polynomial is positive ^ < 0.309. For greater 
value of a this high boundary of ^ tends to constant (see Fig0JV)). Indeed, in the limit a ^ 00 we have 



P 



-5e 



5.2 1. 1 

^ H 

2^ 2^ 21 



(60) 



and it is positive for all ^ < 0.254. We would like to note that for a > 1.136 the polynomial is positive for ^ = 1/6. 
Let us consider now what condition gives the Einstein's equations. The energy corresponding for this configuration 

is 



E = - 



2G 



G\r^{p)dp = - 



2c^ 
G 



1 - 



15 



18 



(61) 
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Equating this energy with zero point energy 



E' 



ircn 




we obtain the relation 



15 -TT- 
18~ 



■2 1 



1 ~ 



a = lpf{am,a,£^). 



To find parameters of stable wormholes of this kind we have to consider this equation at the minimum of function 
f{am, a, ^). Because the function /(am, a, ^) at the minimum is positive, we conclude that the stable wormhole may 
exist for a < 18/(15 — tt^) = 3.5. For a — 3.5 the polynomial is equal to zero for ^ = 0.278. Therefore the stable 
wormholes with this profile of throat may exist only for r/a < 3.5. This upper boundary depends on the model of 
throat. For example, wormhole with profile of throat 



gives another boundary, namely r/a < 36/(7r^ — 6) = 9.3. Specific value of r, a and region of f may be found by 
numerical calculation of the function f{am, a, 



In the paper we analyzed the possibility of existence of the semi-classical wormholes with metric ^ and throat's 
profile given by Eqs.©,®- Our approach consists of considering two heat kernel coefficients B2 and B3. We 
developed a method for calculation the heat kernel coefficient and obtained general expression for arbitrary coefficient 
in background The first seven coefficients in manifest form for arbitrary profile of throat are given by Eqs. 
(|34|1 and (|35|l . The sufficient condition for existence of wormhole is positivity of both B2 and B3. Some additional 
conditions may follow from t ~ t component of the Einstein's equations. 

The common property of both models is that the zero point energy for small size of wormhole is always positive for 
arbitrary constant ^. This statement is opposite to that obtained for zero length throat model in Ref. 11]. The behavior 
of zero point energy for large size of wormhole crucially depends on the non-minimal coupling ^ and parameters of the 
model. We show that the wormholes with the first profile of throat may be a self-consistent solution of semi-classical 
Einstein's equations if the constant of non- minimal connection ^ > 0.266. This type of wormholes is characterized by 
the only parameter a, which is the radius of wormhole's throat. The space outside of throat polynomially tends to 
Minkowkian and there is no way to define the length of throat. We would like to note that the minimal connection 
^ = and conformal connection ^ = 1/6 do not obey this condition. 

The second model of wormhole's throat (|2Jl is characterized by two parameters r and a. The latter is the radius of 
wormhole's throat and first is the length of throat. It is possible to introduce the length of throat because the space 
outside the throat becomes Minkowskian exponentially fast. Suitable illustration for this statement is Fig0IIII). The 
existence of this kind of wormholes crucially depends on the parameter ^ and ratio of the length and radius of throat: 
a = r/a. The general condition for a follows from Einstein's equations, namely a < 3.5. The wormhole with very 
small parameter a may be self-consistent considered by scalar massive field with large value of ^ ~ 1/a. The scalar 
field with conformal connection ^ = 1/6 may self-consistently describe wormholes with r/a G (1.136,3.5). For 1^ = 
we obtain another interval r/a E (1.473,3.5). 

We would like to note that in the limit of zero length of throat a = r/a — s- there is no connection with results of 
our last paper [lli] where we considered wormhole with zero length of throat. The point is that the model considered in 
that paper is singular at the beginning. The scalar curvature is singular at throat and there is a singular surface with 
codimension one. For this case in Ref. j^] the general formulas for heat kernel coefficients were obtained which can 
not be found as limit case of expression for smooth background jj. The reason of this lies in the following. The heat 
kernel coefficients are defined as an expansion of heat kernel over some dimension parameter which must be smaller 
then characteristic scale of background. For smooth background we may make this ratio small by taking appropriate 
value of expansion parameter, but singular background has at the beginning the zero value of background's scale. 
This leads to new form of heat kernel coefficients. Furthermore, in the limit of large box i? — > 00 in this background 
the coefficient -B5/2 survives and it defines the behavior of energy for large size of wormhole. 

Another interesting achievement of the paper is developing of the zeta- function approach j3| . The radial equation 
in this background IjlSI) can not be solved in close form even for simple profile of throat We obtain the general 
formula for asymptotic expansion of solutions H24I) . using which we find the heat kernel coefficients H34|l . (|35|l in 
general form. After renormalization the zero point energy may be expressed in terms of the S matrix of scattering 



r{p) 



p tan — ha 



T 
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problem H44() . H48(l . More precisely, we need only transmission coefficient sn of barrier (|47|l . (|49() . The point is that 
the radial equation for massive scalar field in the background ^ looks like a one-dimensional Schrodinger equation 
l|46|l for particle with potential H47|l . This potential depends on both orbital momentum I of particle and non- minimal 
coupling constant ^ as well as on the radius of throat a of wormhole. 

In first model the domain of ^ for which the energy may possess a minimum is limited from below. The reason of 
this is connected with the fact that the effective mass 

may change its sign for some ^ limited from below. The same situation occurs in the short-throat flat-space wormhole 
where the scalar curvature is negative, too. This is not the case for second model. The scalar curvature in this 
model 

6y2 -f 4a - Sa^ -f (4j/2 + a{-A + ha)) cosh(^) - 2y (-2 -I- 5a) sinh(^) 
r2 sinh* ( f ) (l - a + y coth( I) ) ^ 

may change its sign depending on the parameters of model. For small a it is negative but starting from a = 4/3 the 
domain around y — Q appears, where the curvature becomes positive. This domain becomes larger for larger values of 
a and for a great enough the curvature is, in fact, positive. It is in qualitative agreement with above consideration. 
Indeed for small values of a (negative scalar curvature) we obtained low boundary for ^ and vise versa, for large 
values of a (positive scalar curvature) we obtained upper boundary for ^. 
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APPENDIX A 



In this Appendix we reproduce in manifest form some expressions, which are rather long to reproduce them in the 
text. First of all let us consider the first five terms of uniform expansion 

5+(x) + 5-(-x) = ln(a,) - \ HSl.r,') + / " 1^ ^ + E -'"'^ 

^ k=Q -^-^ I fc=l -^-l J 

fc=0 "^^^ fc=0 

The coefficients with odd powers of v read (we give the integrands, only) 



•^-1 ^ 1/-'- + 3' 



e (-1 + + Ur^) , (-1 + r^) + ^n^r^ + Sr,^ (-1 + 3r^) + 6r,^r^ - + 4r,5r," 



[.-3] : S, 



8^(l + r,2)^/2 



87^(l + f)^/' 
(-l + rf)V4rfc(l + r,'2)r," 



2rl (l - 8rf -|- H{xf - 2rl^ - + 2r^ ((7 - 31r,'^) r," + 2r^ 

))+2r,^((5-8r,'2)r,"+4r,(^)) 



(4) 



-1 - 11?; 



? + 12r'4 - 4r/'2 + 12r/n 



^ ' ^ 'J 1287^(l-f f)"/2 

+ 4ri (-1 + 8rf + 9rf + Sr,"^ + 8r,'rf + ((7 + 43rf ) + 2rf 
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+ 
+ 



2r,^ (-3 + 17rf - 294rf + 40r^' + SSr^'rf ^) + (4 (l5 - 136rf ) r," + 40r^'^) 
Ar^ (-1 - 5rfc'2 + 120r;,'4 + 23rfc"2 + SSr^Vf ^) + 4r/ ((13 - 168rf ) r^" + ISr^^^^) 



(Ala) 



(4) 



r," (^1 + 24r,:2 + Sr,"^ - A^ir^^) + 8r^ (2 (l + 2r^) r," + 7r,^ 
and below are the coefficients with even power of u: 

K] : Eo = ln(a,)-iln(5^r,4) 



(Alb) 



The functions Sp are defined by relation 



{p=Q J-kR j,^Q 

^ — n J —R „ — n 



Below is the list of first four functions Sp with odd indices, (here r = r(p) and z = 1 + k r [p)) 
-1/2 



1 



^ (-1 + r'2 + 2rr") + ^ {1 - 4rr"} 



25 



-rr ■ 



S3 



-5/2 



3r 



e (-1 + r'2 + 2rr"f + {(-1 + r'^) (l + 12r'2) - 2r (-5 + 8/^) r" 



+ z 
+ 



+ z 



+ z 



(2r"2 + rV^^)) + 4rV(4) } + 1 {l + 24/^ - 16r (l + 2/^) r" + 2>2r'' (r"^ + rV^^)) - IGrV^^) } 

e {-19 (-1 + r'^) r'^ - 3r (1 + 5/^) r" + 2r2 (sr"^ + SrV^^)) } 
i {-r'2 (19 + 120/2) + r (3 + 200/^) r" - 2r2 (l5r"2 + 22rV(3)) + 2rV(^)} 

^/2 ^ ^/2 ^ 2rr") + ^ {l014r'4 + 38rV"2 + r'^ (25 - 832rr") + SSrVV^^) } 



175r 



-11/2 



1989rr'- 
32~~ 



(_3/2 ^ ^^//^ 



200 

-13/2 



12155rr"' 
128 



S5 



-7/2^ 



^3 + /2 ^ 2rr")' + |- { (-1 + ^")' (3 + 112r'') + 12r (-1 + r') (1 + r') (4 + 5/^) r" 
- 4r2 (-27r"2 + 48/^"^ - 26rV(3) + 26/^(3)) + 24r=^ (-2/'^ - 2rV"r(3) - r^^) + r'^r^^)) 



64 

8r2 (^-54r"2 + ISGr'V'^ - 69rV(3) + SOrV 



'3^(3) 



8r'Ml2r"-V32rW-') + llr(4) 



- 16r^ (llr(3)2 ^ i5/v(^) + 3rV(^)) + ISrV^*^)} + ^ {l + 112^'^ + 960/^ - 4r (9 + 424/^ + 192/^) - 
+ 16r2 (21r"2 + lOer'V"^ + 28rV(3) + 48r'V(3)) - Mr'' [ir"' + 27rV"r(3) + r'^) + 6/^^) 
+ 32r4 (9r(3)2 + 13r"rW + 4rV(^)) - 32rV('5) } 



[5^2 (_i + /2^2rr") 



X (-29/" + r'2 (29 - 15rr") + 3rr" (-1 + 2rr") + 20rW(3)) + 4 ~ ^'^^ ^'^ ^"^'^^ 2628/^) 

+ r (-15 - 2839/2 + 2592/"^) /' + ^r^ ^128/'^ + 643/^/'^ + 19l/r(^) + 284/^r(^)) 



25 



- 2r^ (l50r"3 + 638r'r"r^^^ + 5r^^^ + 165r'VW) + 8r^ (jr^""^^ + 12r"rW + TrV^^)) } 
'2 (145 + 5256/^ + 20160r'^) + 3r (5 + 1776r'^ + 17952/^) r" 

(^113r"2 + 8440/^"^ + 166r'r^^^ + 3864r'-V(3)^ + ir^ (676r"^ + 2816rV"r(^) + 5r('') + 756r'V('*)) 



+ — i-r 
64 



+ 2 



-11/2^ 

16 



- 8r" (58r(3)2 + gSr'V^^) + 46rV(5)) + 8rV«^) }" 
+ I { (-1 + r-'^) r'^ (175 + 15054r'^) + 2rr''^ (4413 + 1976/^) r" 

- Ur^ (l9r"2 + 813r'2r"2 + 28r'r(3) + 414r'\(3)) + Ar^ (l33r"3 + 638r'r"r(=^) + 221r'V(4)) } 

+ (175 + 30108r'2 + 393408r"') - 8rr'^ (2119 + 88298/^) r" 

64 

+ 4r^ (l33r"'^ + 65998r'V"^ + 196rV(^) + 31264r'V(^)) - 16r^ (628r"^ + 2729rV"r(^) + 791r'VW) 

+ 8r* (69r(3)2 + llO/'r^^) + 54r'r(^)) }] + z' 

+ 2r^ (9r"^ + 5r'A^^^ | + ^ {-r"* (8177 + 332178r'2) + rr'^ (l989 + 396718 

- 4rV2 (21015r"2 + 10168rV(3)) + 2r^ (631r"3 + 2782rV"r(3) + 815/^^'')) | 



r'2 - 9r (1 + 5r'^) r" 



5r'2) r" 



+ z 



-15/2 



13rr' 
128 



'2 r 



12155^r'2 (-1 + r'^ + 2rr") + ^ {r'^ (12155 + 2052348r'^) - 1484372rr'V" 

8 



+ Ar^ (34503r"2 + 16880rV(3)) | 



+ z 



-17/2 



3727125rr"^ 
512 



(-3/2 + rr") 



+ z 



-19/2 



7040125rr 



/6 



1024 



(A2a) 



and here is the list of first three functions Sp with even indices (r = r{R) and z = 1 + k'^r'^{R)) 
So 



-r^z-\ 



82 = -z-'^2r^ 



- (1 - 4rr") + ^ (-1 + + 2rr") 



15r r 



2 ^'2 



Si 



- ^-^r^ [-3r'2 + rr"] - z' 
[4^2 ^ ^/2 ^ 2rr")' + ^ {(-1 + (l + 6r'2) - 2r (-4 + 5/^) r" - 2r^ (4r"^ + rV^^)) + 2r^r^^^ } 
+ ^ {1 + 12r'2 - 4r (3 + 4r'^) r" + 8r^ (sr"^ + 2r'r^^^) - SrV^ } 
- z-^3r^ (-11 (-1 + r'2) r'2 - 2r (l + 5/^) r" + (^r"^ + 5r'r(3))) 
+ i |_ (1 + 4r'2) (-1 + 15r'2) + 2r (l + 53/^) r" +r^ (-17/'^ - 22r'r(3)) +rV*)} 



30Cr'2 (-1 + r'2 + 2rr") + i {3r'2 (5 + 172/^) - 432rr'V" + Ar'' (5/'^ + 7r'r^^^) } 



565r r 



2^/2 



(-3/2 +rr") 





■1395^2^/4- 
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(A2b) 



Below are the functions Wa, cUb and Wc with definitions of the corresponding integrals. 



--a-F32a-' 
o 



dvv 



e^'^'^^ + l 



■In 



-16(7 + 21n2)C+-(3 + 27) 



l + \/l-z^2_ 
- 87r2aln(2a)(l - tanh(7ra))(l - 2^ 



32T:^a(yi + -T:^a[Vi + 6V2 - 47raV3 + 27raV4 - STr^a^^s], 



Here we introduced five functions for cj^ 



(A3) 



ira 1^ ln(2a:^) ^ vdv 
2 Jo cosh^ {-Kau) Jo 



o27raiy 



+ 1 



In 



+ 



1 + Vi - 1 + 



-f 

Jo 



dv 



In 



cosh^ {-Kav) 

cosh^ [ivav) ' 
^ smh{'Rav)dv 
cosh^ {-Kav) 



1 + 

V 



arccos 



and six for 

Us 
where 



229n 
' 2520 



— wa fi {au) arccos udu, 
Jo 

/ /siaJ^j hTra / 

Ji ^ Jo 



1 + 



p di. \ 1 

2 io 



fr,{av)ydv 



cosh (TTflz/) 

4 



i + \/r 



A 1/2 cosh^(7raz/) 4 Jg 



cosh^ {nai/) 4 Jg cosh^(7raz/) [ + 



2 /-oo 



2 /-l 



8 Ji z/'^cosh (Trai/) 16 Jq cosh (Traz/) 



In 
In 



( - 



z/ 



1 + VT 



3(1 + vr^^)2 



(14595n2(i/) - 43638n3(i/) + 47736n4(i/) - 17680n5(z/)), 

87r 



4-77 

315 

^^(29n2(z/) - 10n3(z/)) - --^(32445n2(i/) - 64113n3(i/) + 55692n4(z/) - 17680n5(i/)), 

_^(16n2(r/) - 503(1/)) + ^(6720n2(i/) - 10773n3(j/) + 7956n4(i/) - 2210n5(!/)), 
2297r 



630 ' 
37677 _ 59487r 

~ 3 315 ' 

,2 54477, 88247r 

= 327r£^ £H , 

^ 3 ^ 315 ' 



and 



n2(!/ 
03(1/ 

n4(!/ 

n5(!y 

5'c(i/ 



n\Sc{v)+^-vSc'{v)), 

^\lsc{v) + ^i^5c'(i/) + ^y''Sc"{v)l 

7r^(25c(^) + |g-5c'(.) + |i-^5c"(^)) + ±-^.^Sc"'{v)), 



cosh (7raz/) 
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